home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / claein.z / claein
Encoding:
Text File  |  2002-10-03  |  4.1 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAEEEEIIIINNNN((((3333SSSS))))                                                          CCCCLLLLAAAAEEEEIIIINNNN((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAEIN - use inverse iteration to find a right or left eigenvector
  10.      corresponding to the eigenvalue W of a complex upper Hessenberg matrix H
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CLAEIN( RIGHTV, NOINIT, N, H, LDH, W, V, B, LDB, RWORK, EPS3,
  14.                         SMLNUM, INFO )
  15.  
  16.          LOGICAL        NOINIT, RIGHTV
  17.  
  18.          INTEGER        INFO, LDB, LDH, N
  19.  
  20.          REAL           EPS3, SMLNUM
  21.  
  22.          COMPLEX        W
  23.  
  24.          REAL           RWORK( * )
  25.  
  26.          COMPLEX        B( LDB, * ), H( LDH, * ), V( * )
  27.  
  28. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  29.      These routines are part of the SCSL Scientific Library and can be loaded
  30.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  31.      directs the linker to use the multi-processor version of the library.
  32.  
  33.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  34.      4 bytes (32 bits). Another version of SCSL is available in which integers
  35.      are 8 bytes (64 bits).  This version allows the user access to larger
  36.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  37.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  38.      only one of the two versions; 4-byte integer and 8-byte integer library
  39.      calls cannot be mixed.
  40.  
  41. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  42.      CLAEIN uses inverse iteration to find a right or left eigenvector
  43.      corresponding to the eigenvalue W of a complex upper Hessenberg matrix H.
  44.  
  45. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  46.      RIGHTV   (input) LOGICAL
  47.               = .TRUE. : compute right eigenvector;
  48.               = .FALSE.: compute left eigenvector.
  49.  
  50.      NOINIT   (input) LOGICAL
  51.               = .TRUE. : no initial vector supplied in V
  52.               = .FALSE.: initial vector supplied in V.
  53.  
  54.      N       (input) INTEGER
  55.              The order of the matrix H.  N >= 0.
  56.  
  57.      H       (input) COMPLEX array, dimension (LDH,N)
  58.              The upper Hessenberg matrix H.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAEEEEIIIINNNN((((3333SSSS))))                                                          CCCCLLLLAAAAEEEEIIIINNNN((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      LDH     (input) INTEGER
  75.              The leading dimension of the array H.  LDH >= max(1,N).
  76.  
  77.      W       (input) COMPLEX
  78.              The eigenvalue of H whose corresponding right or left eigenvector
  79.              is to be computed.
  80.  
  81.      V       (input/output) COMPLEX array, dimension (N)
  82.              On entry, if NOINIT = .FALSE., V must contain a starting vector
  83.              for inverse iteration; otherwise V need not be set.  On exit, V
  84.              contains the computed eigenvector, normalized so that the
  85.              component of largest magnitude has magnitude 1; here the
  86.              magnitude of a complex number (x,y) is taken to be |x| + |y|.
  87.  
  88.      B       (workspace) COMPLEX array, dimension (LDB,N)
  89.  
  90.      LDB     (input) INTEGER
  91.              The leading dimension of the array B.  LDB >= max(1,N).
  92.  
  93.      RWORK   (workspace) REAL array, dimension (N)
  94.  
  95.      EPS3    (input) REAL
  96.              A small machine-dependent value which is used to perturb close
  97.              eigenvalues, and to replace zero pivots.
  98.  
  99.      SMLNUM  (input) REAL
  100.              A machine-dependent value close to the underflow threshold.
  101.  
  102.      INFO    (output) INTEGER
  103.              = 0:  successful exit
  104.              = 1:  inverse iteration did not converge; V is set to the last
  105.              iterate.
  106.  
  107. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  108.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  109.  
  110.      This man page is available only online.
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.